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Abstract Chemical Fricke dosimeter in the laboratory can be submitted 
to gamma radiation at low temperatures to study the evolution of oxidation 
reactions induced by radiation, a key process to understand the formation of 
complex molecules. Products generated by the interaction of the different 
elements under radiation can be determined through a mathematical model that 
considers chemical reactions as coupled nonlinear ordinary differential equations 
involving the mass balance of all the species in the reaction. In this paper is 
implemented an alternative way of solving this system of equations, species’ 
concentrations are calculated through an agent-based model implemented in 
Python. The model is a modified version of the prey-predator model where 
each chemical specie involved is considered as an agent that can interact with 
other specie with known reaction rates leading to production (source terms) and 
to destruction (sink terms). Here, the radiation is a factor that affects product 
formation while the bath temperature modifies the reaction speed. This model 
can reproduce experimental concentrations of products and the consumption of 
ferrous ions from a laboratory reaction of irradiation of iron salt solutions at 3 
different temperatures (dry ice, liquid nitrogen, and room temperature). 


Keywords: Chemical reactions; Fricke dosimeter; Prey-predator model; Agent- 
based model. 


1. INTRODUCTION 

In the laboratory is common to use radiation induced chemical oxidation 
reactions to study the formation of biologically relevant molecules in space 
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bodies due to the energy transfer by radiation [1-6]. Chemical analysis, pulse 
radiolysis, and competition kinetic methods have been employed experimentally 
to obtain the large number of rate constants involved in these reactions 
[7-11]. Mathematical models can follow the dynamics of the radiation induced 
chemical reactions using a system of coupled nonlinear ordinary differential 
equations [12-17]. In general, these models cannot be solved exactly because 
the non-linear coupling of the equations that produces a strong dependence 
on initial conditions [13], however there are some numerical trials of solution 
[12,14,17]. There is an effort to make computational simulations to understand 
chemical reactions evolution, using finite-differences [18], Monte-Carlo [19] 
or agent-based models [20,21], Continuing the trials with agent based models, 
here we propose a modified version of [21] considering that each chemical 
species involved is an agent that can interact with other species with known 
reaction rates leading to production and destruction terms as in a prey-predator 
model [22-24]. The radiation is considered as a factor that affects product 
creation while the bath temperature modifies the reaction speed. This model 
can reproduce the formation of products and the consumption of ferrous ions 
in a laboratory experiment. 

2. EXPERIMENTAL SETUP 

Iron salt oxygenated solutions (Fe 0 S0 4 -711,0) in sulfuric acid were irradiated 
at a dose rate of 0.93 Gy/min in a source of Cobalt 60 (following the setup of 
Negron-Mendoza, et al. [2]). Samples were inside Pyrex cells surrounded by 
a Dewar flask filled with different cooling agents to made a thermal bath at 
77 K (liquid Nitrogen), 198 K (dry ice), and 298 K (water at room temperature). 
After irradiation, the samples were melted, and the temperature was adjusted 
to room temperature (298 K). Then, the analysis was performed by UV-Vis 
spectroscopy at 308 nm in a Varian Cary 100 spectrophotometer (California 
USA) Fe 3+ concentrations were determined and compared with the agent- 
based model results. 

3. MATHEMATICAL MODEL FOR OXIDATION REACTIONS 

Chemical oxidation reactions for Fe 2+ in a Fricke dosimeter are well known 
[14,20]: 
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In this case, the last chemical reactions involve 14 species (Fe 2+ , Fe 3+ , H,0, 
H*, H+, H 2 , 0 2 , ‘OH, OH , H 2 0 2 , H 3 0+, HO*, H0 2 , and e aq - (hydrated 
electrons), and seven known reaction rates (r). Standard concentrations for 
the Fricke dosimeter for Fe,SO ( -7H n O were employed, with the oxidation of 
ferrous ions taking place in an acid medium (H,,S0 4 ). Initial concentrations 
were 0.005 M Fe 2+ , 0.012 M 0 2 , 0.400 M H.O 1 , and the rest of the species 
had zero concentration (in agreement with [7]). Reaction rates are: r ~2.0*10 10 
M-’s- 1 , r ~2.3*10 10 M-V 1 , r ~3.0*10 8 M^s" 1 , r 4 ~50 M^s" 1 , r~ r fi ~7.3*10 5 
IVMs -1 , and r 7 ~10 5 M _1 s _1 . 

Kinetics of the radiation induced chemical reactions in a thermal bath 
(eqs. 1) involves the mass balance for each of the 14 species. This balance, 
assuming that is due with the interaction between pairs, can be described by 
coupled nonlinear differential equations considering external source (radiation 
and bath temperature), production (source terms) and destruction (sink terms), 
as proposed on a previous work [15, 19]: 


dN,{t) 

dt 


f +EE'f5* J (<K (') - K WE'2". (')■ 


( 2 ) 


where denotes the molar concentration of the -species at time t,f. is an external 
source, and denotes the reaction rate constant for the two-reactive species 
j, k that produce the specie (i). 

Here, we consider that the external source term f , is the number of 
molecules of the /-species produced in the radiolysis process induced by the 
gamma radiation absorbed in the medium at the thermal bath temperature T. As 
first approximation, we will take it given as a Maxwell-Boltzmann distribution 
for temperature, and as a linear function of the radiation dose intensity I , 
modifying the expression of previous publications [15, 19]: 
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where M is the molecular mass of the species i, M mo is the molecular mass 
of water (18.02 g/mol), N A is the Avogadro number (6.022 x 10 23 mol), is 
the Boltzmann constant (1.38 x 10 -23 m 2 kgs~ 2 Kr‘), and G. a constant, the 
radiochemical yield that gives the number of species i produced when 100 eV 
are absorbed by the system. For the case of Fe 3+ , G. is 15.6 molecules. 

Chemical reactions described by the system of eqs. (1) can be rewritten 
as a system of nonlinear differential equations using eqs. (2,3), one for each 
species involved: 
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This system of coupled nonlinear differential equations can be solved by 
finite difference methods as Sanchez-Mejorada, et al. [14] or Cruz-Castaneda 
et al [15] done. The problem with this solution scheme is that the system has 
a strong dependence on initial conditions and to guaranty convergence some 
extra parameters need to be added to the model. Moreover, there are some 
computational troubles due to the very wide range of the reaction rates that 
differ more than 8 orders of magnitude. 

4. AGENT-BASED MODEL 

An alternative numerical solution for eqs. (4) can be found using an agent- 
based model with sources and sinks that can be interpreted as prey and 
predators. The agents are in a thermal bath at temperature T under a gamma 
radiation immerses in water with free electrons. We made a Python program to 
follow the molar concentrations of each species. The code initializes a region 
of 1,024 x 1,024 pixels, with 40,000 randomly allocated points associated 
with H-.0+, 1,200 points of O,, and 500 of Fe 2+ . In contrast with a previous 
work [20], for every step, the specie of each point of the grid moves with a 
velocity produced by the thermal agitation of the molecules using a Maxwell- 
Boltzmann distribution: 



(5) 


In each grid point is one of the species involved; arrows correspond to velocity 
of the specie given by eq. 5. Two points react in agreement with eqs. (1) when 
they are in the same point. 

After the specie moves, if in the grid point there are two species that can 
react in accordance to eqs. (1), the reaction takes place consuming the “prey” 
elements by the “predators” produced (see Figure 1). If the velocity of a given 
specie takes it out of the border of the network, the specie stays in its actual 
position. The rate of reaction is determined by . After actualizing the entire net, 
the procedure is repeated until the steady-state is reached. In the program, we 
can follow the molar concentration of each species as a function of time (see 
Figure 2). In Figure 3, the molar concentration of Fe 3+ is plotted and compared 
with experimental data. 
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Figure 1: Agent based model for Fricke dosimeter. 



Figure 2:. Molar concentrations of Fe 2+ and Fe 3+ as function of time. 



Figure 3: Molar concentration of Fe 3+ . The points correspond to experimental values 
and the continuous line to the agent-based results, at 77 K (liquid Nitrogen), 198 K 
(dry ice), and 298 K (water at room temperature). 
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5. CONCLUSION 

In this paper we described the kinetics of the oxidation of ferrous ions induced 
by gamma radiation in a Fricke dosimeter in a thermal bath at 77 K (liquid 
Nitrogen), 198 K (dry ice), and 298 K (water at room temperature). Instead of 
using traditional methods to solve the system of coupled non-linear ordinary 
differential equations that describe the mass balance of all the species involved 
(eqs. 4), here we implemented an original agent based model in Python that 
includes in eqs. (4) all of the reaction rates as source and sink terms (like 
in a prey-predator model) with an external term due to the radiation induced 
product formation and the reaction velocity given by a Maxwell-Boltzmann 
distribution. The agreement of the results from our agent-based model with the 
experimental results suggests that our numerical model is reliable for analyzing 
the kinetics of chemical oxidation reactions induced by radiation in a thermal 
bath. This model can be used to understand the kinetics of reactions induced 
by radiation when the temperature of the irradiated medium decreases below 
the freezing point. Further studies and experiments are needed to validate and 
establish the applicability of our model. 

ACKNOWLEDGMENTS 

Partial financial support was given by UN AM under project PAPIIT IN 106215, 
and CONACYT grant 168579. 

REFERENCES 

[1] Negron-Mendoza, A, & Ponnamperuma, C. (1982). Prebiotic formation of 
higher molecular weight compounds from the photolysis of aqueous acetic 
acid. Photochemistry and Photobiology, 36(5), 595-597, doi: 10.111 l/j.1751- 
1097.1982,tb04421 .x. 

[2] Draganic,Z.D., Vujosevic, S. I., Negron-Mendoza, A., Azamar, J.A., & Draganic, 
I.G. (1985). Radiation chemistry of a multicomponent aqueous system relevant 
to chemistry of cometary nuclei. Journal of molecular evolution, 22(2), 175-187, 
doi: 10.1007/BF02101695. 

[3] Negron-Mendoza, A., Albarran, G., Ramos, S., & Chacon, E. (1995). Some 
aspects of laboratory cometary models. Journal of Biological Physics, 20(1), 
71-76, doi: 10.1007/BF00700422. 

[4] Ehrenfreund, P., & Charnley, S. B. (2000). Organic molecules in the interstellar 
medium, comets, and meteorites: A Voyage from Dark Clouds to the Early 
Earth. Annual Review of Astronomy and Astrophysics, 38, 427-483, doi: 10.1146/ 
annurev .astro.38.1.427. 


Radiation Induced 
Oxidation 
Reactions of 
Ferrous Ions: An 
Agent-based Model 


21 





Rivera, A.L. 
Ramos-Bemal, A.S. 
Negron-Mendoza, A. 


[5] Colin-Garcfa, M., Negron-Mendoza, A., & Ramos-Bernal, S. (2009). Organic 
material formed from gamma irradiation of frozen HCN as a cometary ice 
analog: implication to the origin of life. International Journal of Astrobiology, 9, 
279-288. 

[6] Taquet,V.,Wirstrom,E.S.,&Charnley,S.B. (2016). Formation and recondensation 
of complex organic molecules during protostellar luminosity outbursts. The 
Astrophysical Journal, 821(1), 46, doi: 10.3847/0004-637X/821/1/46. 

[7] O’Donnell, J.H.O., & Sangster, D.F. (1970). Principles of Radiation Chemistry. 
Elsevier, New York. 

[8] Schmidt, K.H. (1972). Electrical conductivity techniques for studying the kinetics 
of radiation-induced chemical reactions in aqueous solutions. International 
Journal for Radiation Physics and Chemistry, 4(4), 439-468, doi: 10.1016/0020- 
7055(72)90020-4. 

[9] Castillo, S. Negron-Mendoza, A., Draganic, Z.D., & Draganic, I.G. (1985). The 
radiolysis of aqueous solutions of malic acid. Radiation Physics and Chemistry, 
26(4), 437-443, doi: 10.1016/0146-5724(85)90232-8. 

[10] Draganic, Z.D., Negron-Mendoza, A., Navarro-Gonzalez, R., & Vujosevic, S.I. 
(1987). The presence of polymeric material in radiolyzed aqueous solution of 
ammonium bicarbonate. International Journal of Radiation Applications and 
Instrumentation. Part C. Radiation Physics and Chemistry, 30(4), 229-231, doi: 
10.1016/1359-0197(87)90126-3. 

[11] Bjergbakke, E., Draganic, Z. D., Sehested, K., & Draganic, I. G. (1989). 
Radiolytic products in waters. Radiochimica Acta, 48(1), 73-78, doi: 10.1524/ 
ract.1989.48.12.73. 

[12] Keene, J.P. (1960). Kinetics of radiation-induced chemical reactions. Nature 
188(4753), 843-844, doi: 10.1038/188843b0. 

[13] Malchow, H., & Petrovskii, S.V. (2002). Dynamical stabilization of an unstable 
equilibrium in chemical and biological systems. Mathematical and Computer 
Modelling, 36(3), 307-319. 

[14] Sanchez-Mejorada, G., Frias, D., Negron-Mendoza, A., & Ramos-Bernal, S. 
(2008). A comparison between experimental results and a mathematical model 
of the oxidation reactions induced by radiation of ferrous ions. Radiation 
Measurements , 43(2), 287-290, doi: 10.1016/j.radmeas.2007.11.038. 

[15] Cruz-Castaneda, J., Negron-Mendoza, A., Frias, D., Colm-Garcfa, M., Heredia, 
A., et al. (2015). Chemical evolution studies: The radiolysis and thermal 
decomposition of malonic acid. Journal of Radioanalytical and Nuclear 
Chemistry, 304(1), 219-225, doi: 10.1007/s 10967-014-3711-z. 

[16] Semenov, S.N., Kraft, L.J., Ainla, A., Zhao, M., Baghbanzadeh, M., Campbell, 
V.E., et al. (2016). Autocatalytic, bistable, oscillatory networks of biologically 
relevant organic reactions. Nature, 537(7622), 656-660, doi: 10.1038/ 
naturel9776. 


22 





[17] Morii, Y., Terashima, H., Koshi, M., Shimizu, T., & Shima, E. (2016). ERENA: 
A fast and robust Jacobian-free integration method for ordinary differential 
equations of chemical kinetics. Journal of Computational Physics, 322,547-558, 
doi: 10.1016/j.jcp.2016.06.022. 

[18] Garvie, M.R. (2007). Finite-Difference schemes for reaction-diffusion equations 
modeling Predator-Prey interactions in MATLAB. Bulletin of mathematical 
biology 69(3), 931-956, doi: 10.1007/s 11538-006-9062-3. 

[19] Zhdanov, V.P. (2002). Monte Carlo simulations of oscillations, chaos and pattern 
formation in heterogeneous catalytic reactions. Surface Science Reports, 45(7), 
231-326, doi: 10.1016/S0167-5729(01)00023-1. 

[20] Rivera, A.L., Negron-Mendoza, A., & Ramos-Bernal, S. (2016). “Agent-based 
model of oxidation reactions of ferrous ions”. Journal of Nuclear Physics, 
Material Sciences, Radiation and Applications 4(1), 149-157, doi: 10.15415/ 
jnp.2016.41015. 

[21] Ausloos, M., Dawid, H., & Merlone, U. (2015). Spatial interactions in agent- 
based modeling. In Complexity and Geographical Economics (pp. 353-377). 
Springer International Publishing. 

[22] Berryman, A.A. (1992). The origin and evolution of predator-prey theory. 
Ecology 73(5), 1530-1535, doi: 10.2307/1940005. 

[23] Pascual, M. (1993). Diffusion-induced chaos in a spatial predator-prey system. 
Proceedings of the Royal Society of London B, 251(1330), 1-7, doi: 10.1098/ 
rspb. 1993.0001. 

[24] Pal, D., & Mahapatra, G.S. (2016). Dynamic behavior of a predator-prey 
system of combined harvesting with interval-valued rate parameters. Nonlinear 
Dynamics, 83, 2113-2123, doi: 10.1007/sll071-015-2469-3. 


Radiation Induced 
Oxidation 
Reactions of 
Ferrous Ions: An 
Agent-based Model 


23 





